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The exchange-correlation potentials stemming from the local-density approximation and several 
generalized-gradient approximations are known to have incorrect asymptotic decay. This failure is 
independent of the dimensionality, but so far the problem has been corrected - within the mentioned 
approximations - only in three dimensions. Here we provide a cured exchange-correlation potential 
in two dimensions, where the applications have a continuously increasing range in, e.g., semicon- 
ductor physics. The given potential is a generalized-gradient approximation, which is as easy to 
apply as the local-density approximation. We demonstrate that the corrected potential agrees very 
well with the analytic result of a two-electron quantum dot in the asymptotic regime, and yields 
plausible exchange-correlation potentials for larger two-dimensional systems. 
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I. INTRODUCTION 

Developments in density- functional theory-^ (DFT), 
particularly in the derivation of functionals for exchange- 
correlation (xc) energies and potentials^ can often be 
characterized by a compromise between accuracy and ef- 
ficiency. In quantum chemistry and condensed matter 
physics, applications of generalized-gradient approxima- 
tions (GGAs) are still among the most popular choices 
to evaluate the xc terms in the Kohn-Sham (KS) scheme. 
Even the local-density approximation (LDA) has well 
preserved its popularity due to the several beneficial 
properties, e.g., the computational efficiency, the natu- 
ral compatibility between the exchange and correlation, 
satisfaction of the sum rule of the xc hole etci 

One of the most prominent problems in the LDA and 
in many GGAs is the wrong asymptotic behavior of the 
xc potential. Essentially, this is due to the fact that these 
approximations depend in a oversimplified fashion on the 
(exponentially decaying) electron density, leading to the 
wrong exponential decay of the approximate xc poten- 
tial. A simple and efficient cure for this problem was put 
forward by van Leeuwen and Baerends/^ They provided 
a GGA for the xc potential with the proper asymptotic 
behavior, commonly known as the LB94 functional (or 
potential). The LB94 potential is very accurate for, e.g., 
ionization energies and electron affinities, and it has re- 
ceived considerable popularity in quantum chemistry. 

The three-dimensional (3D) LDA and GGAs are not 
applicable for systems having suppressed dimensional- 
ity^— such as the well-known two-dimensional (2D) ex- 
amples, e.g., semiconductor quantum dots (QDs), rings, 
slabs, and quantum Hall systems. Instead, the 2D 
LDA££ and - to the best of our knowledge - the only 
2D GGA2. are suitable, in principle, to deal with those 
systems. However, these functionals suffer from the same 
deficiencies as their 3D counterparts, most importantly, 
the wrong asymptotic behavior. Motivated by this fact, 



we derive here a xc potential which corrects the asymp- 
totic decay and remains computationally efficient due to 
its simple GGA form. In the derivation we essentially 
follow the steps of Ref. 3, Finally we demonstrate the 
usefulness of the model potential by considering a few 
typical 2D applications: two semiconductor QDs and two 
quantum rings. 



II. CORRECTED LONG-RANGE BEHAVIOR 

Following Ref. [1 we start by expressing the total xc 
potential in the KS scheme as 

v xc = ^ A +p^F(x) (1) 

where 



In Eq. (|T|), F(x) together with x can be considered as a 
quantity evaluating the relevance of the density gradient 
as well as its asymptotic behavior. We remind that x 
has been introduced already for the 2D GGA in Ref. 0. 
We observe that the term p 1 / 2 F(x) scales linearly under 
homogeneously scaling. As it is well known, the linear 
scaling applies to v x exactly, but this may not be the 
case for v c . Nevertheless, in the spirit of LB94 we treat 
v c together with v x . This may be justified considering 
that for r — > 00 the leading contribution to the asymp- 
totics of v xc = v x + v c is due to v x ; in other words, v c 
decays faster 4^ It should be noted, however, that the 
the present correction has effect also elsewhere, i.e., not 
only in the asymptotic limit. 

For the function F(x) in Eq. ([1]) we require F(0) = 0, 
so that the total v xc will be given by the LDA whenever 
Vp = 0. In this way, the result for the homogeneous 
electron gas is exactly recovered. This property is useful 
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also when considering inhomogeneous systems, in view of 
the success of the LDA. 

To the leading order, the model should fulfil l 11 ' 12 
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(3) 



in contrast to the exponential decay of the LDA poten- 
tial^ This requirement translates into the following con- 
dition 



F(x) 



1 



-1/2 



(4) 



We remind that far away from the central region of QDs 14 
p~ e - ar \ a>Q (5) 
and thus, using Eq. (0 in Eq. ((2} we find 

x ~ 2arp- 1/2 . (6) 
Equations (j4j and ((6|) lead to 



F(x) 
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Moreover, since 



In*- — , 
we may re-write Eq. ([7]) as follows: 

4 in a; 

The desired behavior is given by a function 

-r-2 

F(x) = 



1 + 4/3xsinh (x) 



where 



sinh 1 (x) = ln(x + yT-j- x 2 ) 



(7) 
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The parameter /3 can be determined, for example, by 
fitting the model potential against the exact xc potential 
in a solvable test system. We consider the effects of /3 
below. 

Apart from the factor of four in the denominator of 
Eq. ©, and the 2D definition of x, the resulting model 
is very similar to the 3D case^ Moreover, as for the 3D 
case, the exact constraints for the exchange-correlation 
potential resulting from the translational and rotational 
invariance are satisfied. 

Finally we point out that it is straightforward to obtain 
a spin- dependent model potential by simply adding a spin 
index to all the quantities but j3. 
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FIG. 1: (Color online) Exchange-correlation potentials of a 
parabolic (uu — 1) two-electron quantum dot calculated from 
the analytic density. The result from the present model (solid 
line) has been calculated with /3 = 0.01 in Eq. (flOl . The 
chosen value is kept fixed in the following examples. The inset 
shows the relative deviation of the integrated model potential 
from the exact result as a function of /3 (note the logarithmic 
scale). 



III. FURTHER ANALYSIS AND 
APPLICATIONS 

Next we test our model for the xc potential in a 2D 
two-electron QD defined by a radial external potential 
^cxtM = u) 2 r 2 /2 with ui = 1 (see Ref. [Til for information 
about the units). The ground-state (singlet) density has 
an exact analytic expression ) 20 ' 21 



p(r) 



■{e r2 (l+r 2 /2) + ^e 



r 2 /2 



7r(V27T + 3) 

I (r 2 /2) + r 2 h{r 2 /2) + r 2 / 1 (r 2 /2)] }, (12) 



where Iq(x) is the zeroth-order modified Bessel function 
of the first kind. The exact KS orbital is given simply by 
<ys(r) = y/ p(r)/2, and the KS equation can be inverted to 
obtain the exact KS potential from 



^KS(0 = £KS + 



V 2 tp(r) 
Mr) ' 



(13) 



where cks = 2 is the KS eigenvalue . 10 ' 20 The exact xc 
potential is then given by 



vlT t (r)=v KS (r) 



t(r) - v H (r), 



(14) 



where vh is the Hartree potential, which is in fact the 
only term having no analytic expression for the given 
analytic density. However, it can be solved numerically 
with a good precision. 

In Fig. [T]we compare the exact xc potential with the 
LDA potential and with our model, i.e., Eqs. (fTJ) and (fTU|) 
with (3 = 0.01. Whereas the LDA potential decays expo- 
nentially, the model potential has the correct asymptotic 



3 



o 




2 4 6 8 10 
r (a.u) 



FIG. 2: (Color online) Exchange-correlation potentials of a 
42-electron parabolic quantum dot calculated with the LDA 
(dashed line) and the present model (solid line). The dotted 
line shows the desired asymptotic behavior (— 1/r). 
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FIG. 3: (Color online) Exchange-correlation potentials of a 
two-electron parabolic (cj = 0.5, ro = 3) quantum ring calcu- 
lated with the LDA (dashed line) and the present model (solid 
line). The dotted line shows the desired asymptotic behavior 
(-1/r). 



behavior in an excellent agreement with the exact result. 
At r = the density gradient is zero and hence the model 
agrees with the LDA as required. 

To analyze the effect of the parameter (3 in Eq. (fTO)) 
on the xc potential, we consider the relative difference in 
the integrated xc potential weighted by the density, i.e., 



AI= fcPrp(r)\v%r t (r) 



Ml 



Jd 2 r p(r)v^(r) 



(15) 



This quantity is plotted as a function of (3 in the inset of 
Fig. Q] We find that the deviation is overall very insensi- 
tive to (3 (note the logarithmic x-axis). We have chosen 
the value (3 = 0.01 (dashed line) that yields a "reason- 
able" shape for the xc potential and a relatively small 
deviation in AI for this test case. The chosen value is 
kept fixed in the following examples. 

Figure [5] shows the xc potentials for a significantly 
larger parabolic QD with 42 electrons. The confine- 
ment strength is now u) = 0.5 (Ref. fl9l ). Here we have 
performed a self-consistent LDA calculation using the 
octopus code22 and use the LDA density as an input 
for the model potential. Again the correct asymptotic 
limit seems to be recovered by the model potential, al- 
though the tail cannot be followed as far as in the exact 
two-electron case in Fig. [T] This is due to the numer- 
ical instability in the implementation of the LDA cor- 
relation in the regime where the density is very small 
(< 10~ 16 ). Nevertheless, Fig. [5] confirms the improve- 
ment of the model potential over the LDA for a large 
QD system. 

Finally, we consider quantum rings defined by an exter- 
nal potential v eK t(r) = w 2 (r — ro) 2 /2, where ui = 0.5 and 
the ring radius ro = 3. The obtained potentials for two- 
electron and ten-electron rings are shown in Figs. [3] and 
21 respectively. Similarly to the previous test cases, the 
correct asymptotic behavior (—1/r) is obtained with the 
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FIG. 4: (Color online) As in Fig. [3] but for ten electrons. 



model potential for both quantum rings. At the ring ra- 
dius where the density gradient is zero, the model agrees 
with the LDA in a smooth fashion. 

It is noteworthy that the present model for the xc po- 
tential sets specific constraints on the external potential. 
In particular, infinite potentials (at finite positions) - 
leading to situations where both the density and its gra- 
dient are zero - may result in divergence of F(x). A typ- 
ical example would be a quantum ring studied in Ref. Il2l . 
In this respect, the present potential needs to be applied 
with special care, or alternative models should be consid- 
ered, when studying 2D structures represented by such 
external potentials. Nevertheless, the presented scheme 
opens up a path to design of a more general GGA ap- 
proach. 
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IV. SUMMARY 

In summary, we have derived a generalized-gradient 
approximation for the exchange-correlation potential 
with a proper asymptotic decay for two-dimensional 
nanostructures. The expression has a proper asymp- 
totic decay and thus corrects one of the most signifi- 
cant errors in the commonly used two-dimensional local- 
density approximation. The validity of the model poten- 
tial is confirmed in quantum dots and rings with differ- 
ent number of electrons. We believe that the presented 
approach gives promising prospects in two-dimensional 
many-body physics, first and foremost for systems where 



the correct decay of the exchange-correlation potential is 
of great importance, for example, quantum-dot molecules 
and chains, coupled quantum rings, and periodic electron 
lattices such as the artificial graphene^ 
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